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1 Introduction 



The importance of zeta-function regularisation for the definition of functional determinants Q 
is, without discussion, a powerful tool to deal with the ambiguities (ultraviolet divergences) 
present within the one-loop or external field approximation, in relativistic quantum field theory 
(see the seminal papers |2],||] and for recent reviews @,||]). Let us remember that the Euclidean 
partition related to a quantum scalar field can be formally written as 

Z = (detLc)- 1 / 2 , (f.f) 

with Ly) an elliptic differential operator. The latter quantity is ill defined and, with regard 
to this, we briefly recall how zeta-function regularisation works: it gives a precise meaning, in 
the sense of analytic continuation, to the determinant of a differential operator which, as the 
product of its eigenvalues, is formally divergent. When the (.D-dimensional) manifold is smooth 
and compact, the spectrum is discrete and one has, for Re s > D/l, I being the order of the 
differential operator Lp, the definition of the zeta-function 

c(s\l d ) = j2k s , 

% 

where \ are the eigenvalues of the elliptic operator Lq. Making use of the relationship between 
the zeta-function and the heat-kernel trace via the the Mellin transform, when Re s > D/l, one 
can write 

1 r°° 

C(s\L D ) = TrL D s = — f- 1 K(t\L D ) dt , (f .2) 

1 (s) Jo 

where K(t\Lo) = Tr exp(— tL^) is the trace of the heat operator. The previous relations are 
valid also in the presence of zero modes, with the replacement K(t\Lo) — ► K(t\Lu) — Po, Po 
being the projector onto the set of zero modes. 

A well known heat-kernel expansion argument leads to the meromorphic structure of (,(s\Ld). 
It is found that the analytically continued zeta-function is regular at s = and thus its derivative 
in that point is well defined. As a consequence, the one- loop Euclidean partition function, 
regularised by means of the zeta-function, reads || 

InZ = -Ilndet(^ 2 L D ) = ^'{<d\l 2 L D ) = \c,'(Q\L D ) - ic(0|^)ln€ 2 , 

where £' is the derivative with respect to s and £ is a renormalization scale. 

However, in general things are not so simple. Sometimes it happens that one has to deal 
with the product of two or more differential operators and one is directly confronted with the 
validity (or not) of the multiplicative property: 

lndet(AB) = bidets + In detB , (1.3) 

which is of course known to hold for non singular matrices in finite dimensional vector spaces. 

A first elementary example is the following. Consider a free vector-valued scalar field (j>i in 
R 4 , with a broken O(N) symmetry — owing to the mass terms mf. The Euclidean action is 

S = f dx 4 ^ [(-A + m?) , (1.4) 
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and the related Euclidean operator reads 

= (-A + mf) Sij , (1.5) 



in which A is the Laplace operator. Thus, one is actually dealing with a matrix-valued elliptic 
differential operator. In this case, the Euclidean partition function is given by 



In Z = — In det 
2 



1 _ , 
- In det 
2 



£ 2 (-A + m 2 )..l 2 (-A + m 2 N ) . (1.6) 



Another, less trivial, example concerns the finite temperature effects for a gas of free rela- 
tivistic charged bosons. With regard to this case, it is possible to show that the logarithm of 
the grand canonical partition function, choosing as parametrization of the charged boson field 
the two real scalar fields fa, can be expressed by || 

In Z 0tlt = ~ In det (j^L+L.) , (1.7) 

with 



L± = -d 2 - A + m 2 + [i 2 ± 2fiy/-A + m 2 , (1.8) 

r being the imaginary time compactified with period (3, the inverse of equilibrium temperature, 
and \x the chemical potential. Note however that, if one chooses as parametrization of the 
charged boson field the complex scalar fields 4> and <p*, one has 

hiZfo = -i In det (t^K+xJ) , (1.9) 

with 

K± = -d 2 - A + m 2 + fi 2 ± 2fid T . (1.10) 

In a recent work, Stuart Dowkerjfjj has shed doubt on the recipe used in computing the 
partition functions, Eqs. (|l.6| ) and (|l.7| ). In our opinion, the above prescription is correct and 
we refer to Ref. ||, where this issue has been discussed in detail. 

As further examples, we simply recall that in the evaluation of the one-loop Vilkovisky- 
DeWitt effective action (see, for example, Refs. 0- [12| for details) and in some GUT-like models 



(see, for example, [IH]), one has to deal with matrices of higher order differential operators, which 
give rise to products of the same, when one has to compute functional determinants. 

Within the physical literature, in all the examples we have recalled, and in many other that 
we do not mention here, the way one usually proceeds is by formally assuming the validity 
of the multiplicative relation, indiscriminately. Needless to say, this may be dangerous. One 
has to use always some regularisation procedure and it turns out, in fact, that the regularised 
determinants do not satisfy in general the multiplicative property, Eq. (|b^). Even when one is 



dealing with commuting operators, there exists the so-called multiplicative anomaly [14, 1^]. In 
terms of F(A,B) = det(AB) / (det A det B) pfl, it is simply defined as 

a(A,B) = \aF(A,B) = In det (AB) -lndet(^) - lndet(B) . (1.11) 

It should be noted by passing that the non vanishing of the multiplicative anomaly implies 
that the relation 

In det A = Tr In A (1.12) 
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does not hold, in general, for elliptic operators. The formal use of the above operator identity 
is not justified, since the multiplicative anomaly may be present. In fact, if one assumes that 
Tr In A is a linear functional, in the simplest case of [A, B] = 0, one has 

Tr hiAB= Tr (ln^ + ln.B) = Tr ]nA + Tr In 5. (1.13) 



Thus, if Eq. (P3 ) holds, one arrives at a contradiction with Eq. ( |1.11[) , as soon as the multi 



plicative anomaly is not vanishing. With regard to this issue, taking for granted Eq. ( |1.12| ), one 
might start with the definition 

Z = exp(-^ Tr lnL D ) , (1.14) 



instead of Eq. ( |1.1|) and then make use of some regularisation while preserving the linear property 
of the trace. In this case, in Ref. [|r| the absence of the multiplicative anomaly has been claimed 
and doubt has been shed on the use of zeta-function regularisation. In our opinion, this is not 
justified (see Ref. JTtJ). 

The presence of the multiplicative anomaly in the commuting case might turn out as a 



surprise. However, in Ref. [18| it has been shown that its existence is actually unavoidable in 
order to have independence of the partition function on the choice of bosonic degrees of freedom 
(see, also ||). In some cases, the multiplicative anomaly gives a contribution to the effective 
action which can be readily absorbed into the renormalization procedure, but in some systems it 
may certainly produce physical consequences, which have to be carefully analysed. A non trivial 
example has been discussed in Ref. [|18| , where after a renormalization of the charge operator, it 
has been shown that, starting from the natural definition Eq. fll.lD , the multiplicative anomaly 
gives rise to a new contribution — overlooked in previous treatments — to the high temperature 
expansion of the free energy of the relativistic boson gas in the symmetric unbroken phase. 

We will show, and illustrate with the help of several examples, that this multiplicative 
anomaly appears already in very simple situations (one-dimensional, first order differential op- 
erators differing in a constant term). And that it can be most conveniently expressed by means 
of the non-commutative residue associated with a classical pseudo-differential operator, known 
as the Wodzicki residue p9| . In fact, the purpose of the present paper is to continue the analysis 



of the emergence of the multiplicative anomaly, that was started in Refs. |[|,20] — where the 
commutative case was discussed — and to extend it to the more general non-commutative case, 
working out explicitly the multiplicative anomaly formula for a large class of elliptic operators 
— and to investigate some very basic examples with care. In particular, in some of these ex- 
amples we will deal with non compact manifolds, where, strictly speaking, the Wodzicki theory 
has not yet been developed. Nevertheless, we will see in such examples that the formula for the 
multiplicative anomaly turns out to be valid too. 

The content of the paper is as follows. In Sect. |2| a perturbative derivation of the multiplica- 
tive anomaly for a particular class of differential operators in arbitrary dimensions is presented, 
while the general formula in the case of lower dimensions D < 4 is given in ^4]. In Sect. || 
we shall recall the definition of the Wodzicki residue, together with some related results, which 
will be used in Sect. |] in order to render explicit its relation with the multiplicative anomaly 
formula. Then, we show in |4.1| that such a general formula can be notably simplified in lower 
dimensions, confirming completely the result obtained by using perturbation theory. In Sect. ||, 
we show how the multiplicative anomaly might be used in order to compute the heat-kernel coef- 
ficients. Finally, in Sect. || we discuss some basic physical examples, formulated in non-compact 
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manifolds, but for which the spectrum of the Hamiltonian operator is exactly known, and we 
compute the anomaly explicitly. The paper ends with some conclusions in Sect. [7|. 



2 Perturbative derivation of the the multiplicative anomaly 

In this section we will consider a particular case involving two self-adjoint elliptic invertible 
operators H, Hy = H + V, and the related product 

A = H(H + V). (2.1) 

By definition, the corresponding multiplicative anomaly reads 

a(H, H v ) = ln det A — In det iT — In det Hy (2.2) 

where the functional determinants are evaluated by using ^-function regularisation. It is conve- 
nient to introduce the quantity 

A(s) = C(s\H) + ((s\Hv) - ((s\A) , (2.3) 

where ((z\L) is the zeta- function associated with the elliptic operator L. Thus 

a(H, Hy) = A'(0) = lim -f [((s\H) + ((s\H v ) - ((s\A)] . (2.4) 
s— >o ds 

Now we also suppose V to be a small perturbation potential. We indicate by G = H^ 1 
the inverse operator of H and by Gy the inverse of Hy. One has the well known operatorial 
equation Gy = G — GVGy , whose solution is formally given by 

oo 

Gy = G J2(- VG T ■ ( 2 - 5 ) 
n=0 

To begin with, let us consider the special, but physically important case of a constant V. 
Then H and Hy are commuting operators, and we may compute the complex power by means 
of the binomial expansion, 



G s v = G s 



1 + E Us)V n G r - 



Jn{S) = itTT \ • (2.6) 

nil (s) 



71=1 

Using the definition of the zeta-function and the properties of the trace, we obtain 



C(s\H v ) = Tr G s v = C(s\H) + £ f n (s) Tr (V n H~ s ~ n ) (2.7) 

71=1 

and, in a similar way, 

oo 

C(s\A) = C(2s\H) + Us) Tr {V n H~ 2s ' n ) . (2.8) 
n=l 

Finally 

oo 

A(s) = E fn(s) [C(s + n\H) - C(2s + n\H)\ tr V n , (2.9) 

n=l 
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where tr is the trace on internal indices. Under the assumptions above, all operators here 
involved are invertible. If H has a zero eigenvalue this must be excluded from the definition of 
the zeta functions for the operators H and A, but it gives a contribution to £(s\Hy). Such a 
contribution modifies the anomaly by an additive logarithmic term since, in this case, 

oo 

A(s) = g V~ s + J2 Us) [C(« + n\H) - C(2s + n\H)} tr V n , (2.10) 
n=l 

go being the number of zero modes. 

As is well known []2l|], for elliptic operators the zeta function has only simple poles. Then 
the multiplicative anomaly assumes the form 

a(H, H v ) = A!{Q>) = -g lnV+Y, (~l) n 7 + ^ n) ^s ((s\H)\ s=n tr V n , (2.11) 

n>2 



where 



n— 1 i 

7 + ^) = £t, (2.12) 



k =i k 



the sum in Eq. ( 2.11 ) being extended to all positive integers where the zeta- function has a 
simple pole. Here 7 is the Euler-Mascheroni constant, while ip{x) is the digamma function. This 
expression has limited validity, mainly because it assumes that the operators involved in the 
product commute, but nevertheless, one should note that it holds as soon as one can define the 
zeta-functions involved. That includes, for example, the case when the manifold is non-compact 
or when it is compact and with boundary. 

Furthermore, when one deals with the last case, i.e. a compact D-dimensional manifold 
with boundary, one can certainly proceed, since the meromorphic structure of the zeta function 
corresponding to an elliptic operator H of order h is well known. In fact, the short-t heat-kernel 
asymptotic expansion of Tr e~ tH is known (see, for example, |22j). Making then use of the 
Mellin transform 

1 r°° 

C(s\H) = — - / dtt s ~ l Tr e~ tH , (2.13) 
a standard procedure yields the analytical continuation 

C(s\H) = -i-£ K l_ D + analytic part , (2.14) 

1 \ S ) n s 

where 



h 



K n = J K{x\H)d D x (2.15) 

are the heat-kernel coefficients related to H, which enter in the asymptotic expansion of Tr e~ tH . 
They are, in principle, computable (see, for example, for the case of a Laplace-type operator, 
the recent paper [p3| and references therein). From this we immediately obtain 



Res C(s\H)\ s=n = K ^ U * ] (2.16) 
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and thus, 



a(H,H v ) = -g \nV+ ]T (-1)" - 'ff 1 K D ^ hn (H) tr V n . (2.17) 

n>2,hn<D H ' 

We close this section with some remarks. The multiplicative anomaly formula above is valid 
only for constant V. For a manifold without boundary one has K r = 0, when r is odd. Thus, for 
a second order differential operator (h = 2) in any odd-dimensional compact manifold without 
boundary, the multiplicative anomaly vanishes. It is also vanishing in two dimensions, but it is 
actually present for D = 4, being proportional to the Seeley-DeWitt coefficient Kq. In the case 
of first order differential operators (h = 1), the anomaly is non- vanishing for D = 2, and for the 
physically more interesting case D = 4 one has (assuming zero modes to be absent) 

a(H, H v ) = \k 2 (H) tr V 2 + ^K (H) tr V 4 . (2.18) 

We note the presence of the first non trivial Seeley-DeWitt coefficient K2(H), which for a Dirac- 
like operator depends on the scalar curvature of the manifold one is dealing with and V may be 
a mass difference. This could be interpreted as a (potentially) interesting, finite effect of induced 
gravity -according to Sakharov- caused by the quantum fluctuation of the matter spinor field. 

It is also interesting to note that the linear term in V does never contribute to the multi- 
plicative anomaly. If one invokes dimensional arguments, this fact is not really surprising in 
the case of constant V. However, as we shall see in the following, this result is also true for an 
arbitrary potential. To show this important fact we shall make use of standard perturbation 
theory. 



2.1 Explicit expression of the anomaly from perturbation theory 

As above, we assume H to be an elliptic differential operator of known, non-degenerate spectral 
decomposition {Aj, ¥>i} i€l and V a small (non constant in general) perturbation potential. Let us 
denote by [ii and at the eigenvalues of the operator Hy and A respectively. By using perturbation 
theory up to second order in V, we have 

fii = k + Vu + Y^P^ + OiV 3 ), (2.19) 

a t = Af + A^ + ^ A f 2 ^ +0(y 3 ), (2.20) 

Vn = {<fi,V Vj ) . (2.21) 

Let us assume that Re s is sufficiently large. Then, by definition C{ S \H) = J2i \ rS while, using 
again the binomial expansion 

C(s\H v ) = £/z- s = 5 oF- s + C(s|tf)-s£^A- S_1 (2-22) 

i i 

+£ ( £ +l) ^ vix _ s _ 2 _ s ^ JMl A - S -i + oiV 2 ) , (2.23) 
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C(s\A) = J2 a 7 S = ((MH)-sY; VuK 2s ~ l (2-24) 

i i 

+ f C£+l) ^ ^ ?A _ 2s _ 2 _ g ^ V^ A - 2S -1 + 0{ y 2) . (2>25) 



2 ^ \^.\?-A? 
As a result, the quantity -4(s) can be written as 

-A(s) = 5o^" s + 2((s\H) -((2s\H) + a[*i(2s) - 
s(s + l 



2 

where we have introduced the functions 



[$ 2 (2s) - $ 2 (*)] + *[*i(2«) - + * 2 (2s)] , (2.26) 



*i00 = E ^fV- ' , (2-28) 

*2(.) = E (2 - 29) 

Since we are interested in the derivative at s = 0, let us assume that the behaviour of the 
analytic continuation of these functions, in a neighborhood of the origin, has the form 

*fc(s) = — + A k + 0(s), (2.30) 
s 

**(*) = — + £> fc + 0(s), fc = l,2. (2.31) 
s 

in agreement with the fact that the zeta function has only simple poles. From Eqs. ( 2.26| )- (|2.31[) , 
we easily obtain 

a(H,H v ) = ^- + D 2 , (2.32) 

again up to second order of perturbation theory. As anticipated before, at first order in V 
there is no contribution to the multiplicative anomaly (observe that the potential here is not 
necessarily constant). 

When V is constant, one has 

$k(s) = V k £(s + k\H), * 1 ( s ) = * 2 ( s ) = (2.33) 

and then 

a{H,Hy) = ^ m K D _ 2h {H), (2.34) 

which is the leading term in V in Eq. ( 2.11 ) . To summarize, up to second order of perturbation 
theory and for the special class of product of two operators, the multiplicative anomaly depends 
only on V 2 and not on the derivatives of V. Moreover, for dimensional reasons, in a 4-dimensional 
manifold and for a second order differential operator, the second order approximation gives the 
exact result and thus, in such case one may argue that the form of the anomaly reads 

a(H, H v ) = — ^ J ka(x\H) tr (V 2 ) ofx , (2.35) 

whatever be the potential (e.g. not necessarily constant). We will confirm this result making 
use of a more powerful technique we are going to introduce in the next section. 
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3 The Wodzicki residue 



For the reader's convenience, we will update in this section some basic information concerning 
the Wodzicki residue (see also [14] and the references to Wodzicki quoted therein) that will 
be used in the rest of the paper. 

Let us consider a D-dimensional, smooth (compact) manifold without boundary, Mjj, and 
a (classical) ^DO, Q, acting on sections of vector bundles on Md- To any classical ^DO, 
it corresponds a complete symbol cr(Q) = Q{x,k) = e tkx Qe~ %kx , such that, modulo infinitely 
smoothing operators, one has 

TT^D I dye^ k Q(x,k)f(y). (3.1) 
Jrd {2it) u Jrd 

The complete symbol of Q admits an asymptotic expansion for \k\ — > oo, given by 

oo 

Q(x,k)~J2Qq-j(x,k), (3.2) 

3=0 

where the coefficients fulfill the homogeneity property Q q -j(x,tk) = t q ~ 3 'Q q _j(x,k), for t > 0, 
being Q q (x, k) ^ 0. The number q is called the order of Q. We shall deal mainly with self-adjoint 
operators and thus we will assume that their complex powers, beside the semigroup property, 
also satisfy (A c )~ s = A~ cs , the c and s being arbitrary complex numbers. As an example that 
we will encounter, let us consider In A, where A is an elliptic operator of order a. Then 

oo 

cr(lnA) ~ aln|fc| +^A_ j (x,/c) . (3.3) 
3=0 

Observe that In A is not a classical 'I'DO operator, since in general its symbol differs from that 
of a zero order operator in the presence of the In \k\ term. 

In order to introduce the definition of the non-commutative residue of Q, let us consider an 
elliptic operator A, with a > q, and form the family of ^DOs Aq(u) = A + uQ, u being a real 
parameter. The associated zeta- function reads 

C(s\A Q (u)) = Tr (A + uQ)~ s . (3.4) 

The meromorphic structure of the above zeta-function can be obtained from the short-t asymp- 
totics of Tr e ~ A< 5 (u) [M, namely 



Tr e -tA Q (u) „ aj (u)t^- D ^ a + Pk(u)t k Int . (3.5) 

j=l k=l 

Note the presence of logarithmic terms that lead, using the Mellin transform, to double poles in 
the meromorphic expansion of £(s\Aq(u)), i.e. 



T(s) 
1 



(3.6) 
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where J(s,u) is the analytic part. Taking the derivative with respect to u and then the limit 
u — * 0, one gets 



lim 4- Tr ( A + uQY 

u^o du 



-s Tr (QA 



-s-l\ 



1 



r s) 



E 



«'(0) 



By definition, the non-commutative residue of Q is 



res(Q) = Res 



a lim — Tr (A + uQ)~ 



0^(0) , 



(3.7) 



(3.8) 



where Res is the usual Cauchy residue. It is possible to show that res(Q) is independent on 
the elliptic operator A and that it is a trace in the algebra of classical ^DOs (actually, the only 
trace up to multiplicative constants). From the above definition and taking the derivative with 
respect to u at u = of Eq. (p^q), one obtains a possible way to compute the non-commutative 
residue. In fact 



TY [Qe 



-tA 



i=i 



\nt + 0(t Int) 



(3.9) 



and so the non-commutative residue of Q can be read off from the short t asymptotics of the 
quantity Tr (^Qe~ tA ^j, just picking the coefficient associated with Int. When the manifold is 
non-compact, this is one of the methods that we have at hand for evaluating the Wodzicki 
residue, as long as all the traces involved exist. 

For the case of a compact manifold, Wodzicki has obtained a useful local form of the non- 
commutative residue, that is, a density which can be integrated to yield the non-commutative 
residue, namely 



res(Q) 



dx 



M D (2tt) d J W=1 



Q^d(x, k)dk . 



(3.10) 



Here the component of order —D (remember that D is the dimension of the manifold) of the 
complete symbol appears. 

Let us now consider an elliptic (self-adjoint) operator B of order b > q. From Eq. ( |3.7| ) and 
the requirement Tr [Q(B C )~ Z ] = Tr {QB~ CZ ). One has the following 



Lemma 1. In a neighborhood of z = 0, 
res(Q) 7res(Q) 



Tr (QB~ 



+ 



R Q (B)+0(z) 



zb b 

where 7 is the Euler-Mascheroni constant and the quantity Rq{B) satisfies the relation 

7 res(Q) 1 — c 



R Q (B C )-R Q (B) 



c> 



(3.11) 



(3.12) 



b c 1 

The latter equation gives another way to compute the non-commutative residue, namely 

res(Q) = b [ Res Tr (QB~ z )] z=Q . (3.13) 
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Making use of Lemma 1, one also obtains (with c\ > and C2 > 0): 

Tr [Q(A— - B-™ )] = ( J_ _ M + 7 res(Q) (I _ 1 

z \aci 0C2/ \a b 

-R Q (A)+R Q (B) + 0(z). 
In particular, if aci = 6c2 = 1 there is no pole at z = 0, and then 

Tr [Q(,T5-.B-f)] = 7 res(Q) Q - - Rq(A) + Rq{B) + O(z) . 



(3.14) 



(3.15) 



On the other hand, we also have [15| 
Lemma 2 



Tr 



Q{A-a -B-V 



-res 



Q 



InA InB 



a 



+ 0(z) 



(3.16) 



This Lemma is a consequence of Lemma 1. In fact, if C is a suitable elliptic operator, making 
use of the relation C z = I + zlnC + 0(z 2 ), one obtains 



Tr 



Q(A~a - B~b] 



Tr 



Q(A~a - BS)C Z C- 



-z Tr 



+ 0(z 2 ). 



(3.17) 



Now Qi^-^- — is a classical ^DO and we can make use of Lemma 1 to obtain the desired 
result. It is easy to show that Lemmas 1 and 2 allow us to rewrite Eq. ( p. 15 ) as 



res(Q) 

z \a b 



Tr [Q{A~ Z -B- Z ) 
As a result, we have 
Lemma 3 

lim^jzTr [Q(A--B-*)]} 



res 



Q 



InA InB 



a 



+ 0(z) 



(3.18) 



-res 



Q 



\nA \nB 



(3.19) 



4 The multiplicative anomaly formula 



In this section we will present a quick proof of the multiplicative anomaly formula, following the 
derivation in [JD^]. We consider two invertible, elliptic, self-adjoint operators A and B on Mjj of 
orders a and b respectively, and the quantity 



F(A,B) 



det(AB) 



a(A,B) 



(4.1) 



(detA)(det5) 
Moreover, we introduce the family of ^DOs 

A(t) = r/*sf, A(0) = sf, A(1)=A, 

tj = AB~i , ord?7 = 0, ord^(t)=o, ord (A(t)B) = a + b . (4.2) 
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and the function 
F(A(t),B) 



One trivially gets 

F(A(0),B) = l 



deb(A(t)B) 
(detA(t))(detB) ' 

F(A(1)),B) 



det(AB) 



F(A,B). 



(4.3) 



(4.4) 



(det A)(detB) 

As a consequence, one is led to deal with the following expression for the anomaly 

a{A(t),B) = In F{A(t),B) = - lim d s [ Tr (A(t)B)- s - TV A(t)~ s - Tr B~ s ] . (4.5) 

s— >0 

This quantity has the properties: a(A(Q),B) = and a(A(l),B) = a(A,B). 

The next step is to compute the first derivative of a(A(t),B) with respect to t, the result 
being 



d t a(A(t),B) = lirnd s {s TV In 77 [(A(t)5)~ s - Tr A(ty s ]} . 



Furthermore, making use of the Lemma 3, i.e. Eq. ( |3,19| ), with Q = In 77, one obtains 

■ In A(t) ln[A(t)BY 



d t a(A(t),B) = res 



In 77 



a + b 



(4.6) 



(4.7) 



Finally, performing the integration with respect to t from to 1, one gets the Kontsevich-Vishik 
multiplicative anomaly formula |]15|, namely 



a(A,B) = f 
Jo 







/ dt res 


In 7/ ^ 


/o 



'lnA(t) ln[A(t)S]\ 



a + 6 



(4.8) 



The multiplicative anomaly formula, Eq. (4.8), notably simplifies in the special case of com- 
muting operators. In fact one has 



hx[A(t)B] tlnr] InB 



a + b 



a + b 



+ 



\nA{t) _ tin 77 \nB 



and 



\n[A(t)B] hxA{t) _ btlnr] 



a + b 
As a result, 

a(A,B) 



a(a + b) 



2a(a + b) 



res 



(ln(AS-6)) 2 j , 

which can be rewritten as the Wodzicki multiplicative formula [14] 
\\n(A b B- a )) 2 



(4.9) 



(4.10) 



(4.11) 



res 



a(A,B) 



2ab(a + b) 



a(B, A) 



(4.12) 



where the symmetry property in A and B is manifest. 

We would like to end this section with the following simple (but important) remark. The 
notion of non-commutative residue can be introduced only for classical "fDOs, namely the ones 
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whose symbol admits the asymptotics Eq. ( |3.2[ ). For example, in general In A with A elliptic 
and a > 0, is not a classical ^DO, but — is, as well as In 77, since ord(r/) = 0. Thus, all 
the formulae in which the non-commutative residue appears are well defined. 
For more than two operators, the recurrence equation 

a(A,B,C) = a{AB,C)+a{A,B) (4.13) 

holds and thus, the knowledge of the multiplicative formula for the product of two operators is 
sufficient for computing the multiplicative anomaly associated with the product of an arbitrary 
number of operators. 

4.1 The multiplicative anomaly formula in lower dimensions 

In the commutative case, examples of the evaluation of the multiplicative anomaly have been 



already presented (see [18,2(J). For the non commutative one, we have obtained some insights 
by using perturbative arguments. Here we would like to compute the multiplicative anomaly in 
the lower dimensional cases D < 4, by making use of the Kontsevich-Vishik formula Eq. Q4.8| ), 
We shall consider a compact manifold, in order to make use of the Wodzicki local formula for the 
evaluation of the non-commutative residue. With this aim, we start from the following property 
for the product of the two non-commuting self-adjoint elliptic ^DOs that we shall consider, Q 
and P (Baker-Campbell-Haussdorf formula) 

In (PQ) = InP + In Q + [In P, G] + [In Q, [In Q, In P]] + 0(ln 3 Q InP) , (4.14) 

where G is a suitable nested commutator. Making use of this operator identity, we have that 
the zero order ^DO Qt, defined by 

Qt = -t — , 4.15 

a a + b 

may be rewritten as 

Qt = - Hln \ + - r^r[lnB, [In B, In 7?]] + 0(ln 3 B In rj) . (4.16) 

a(a + 0) 126 



Since res (In TjQt) is the integrand in the multiplicative anomaly formula, Eq. (|4.8| ), and the 
non commutative residue is a trace in the algebra of classical ^DOs, one has 

res (In t? [In r/, G]) =0. (4.17) 

As a consequence, the Kontsevich-Vishik formula Eq. ( f4.8| ) yields 



res (in 2 77 ) 1 
<A,B) = 2ab K + ^ - —res (lnr?[lnP, [mi?, In 77]]) 

+0 (res (In 2 77 In 3 B)\ , (4.18) 
in which the term which is non vanishing in the commutative case has been conveniently isolated. 
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Now we recall that for a zero order ^DO, rj, one has a(F(j])) = F{a{rf)) (for F{.) analytic) 
and, moreover, that the symbol a of the commutator of two ^DOs reads 



dx^ 



' 9 



dk 



+higher order derivatives 



(4.19) 



If we assume, as in Sect. ||, that A and B have the same principal and sub-leading symbols 
except for the homogeneous one of degree zero, namely a = b and cr(A) — cr(B) = V(x), then 



a(hxr,) = 0(l/\k\ a ), 



d 



a(lnr,) = 0(l/\k\ a ), 



8 



dxt 1 

and it follows that 



<r(ln5) = 0{\kf 



d 
8k» 



a(lnB) = 0(l/\k\), 



(4.20) 
(4.21) 



a(]nr][hiB, [lnB,lnry]]) = 0(l/\k\ 2a+2 ) . (4.22) 
This means that if the dimension of the compact manifold D satisfies D < 2a + 2, only the first 



term on the right-hand side of Eq. ( 4.18 ) gives a non vanishing contribution to the anomaly. 
Such a condition is satisfied in the particular but important case of second order differential 
operators in four dimensions. On the other hand, for first order differential operators a = b = 1 
and D = 4, also the second term of Eq. ( 4.18| ) gives a contribution. This exact analysis confirms 
the perturbative results of Sect. |2[ 



5 Heat kernel coefficients from the multiplicative anomaly 

As we we have shown in Section |2l Eq. ( 2.17 ), in the case of constant V the multiplicative 
anomaly a(H,Hy) can be directly related to the heat kernel coefficients K n {H), but it can 
also be computed by means of Wodzicki formula, Eq. ( 4.12 ). This means that one could use 
the Wodzicki residue in order to compute heat-kernel coefficients (see Ref. [14,^]). Now we 
illustrate the method for an invertible differential operator H, of order h, on a manifold without 
boundary. 

For this case, Eq. ( 4.12] ) reduces to 



1 r _ 1 -i2 1 

ci£)(Hv,H) = — res IniHyH ) = — res 
4aj l J 4n 

If V is a constant, from Eq. ( 2.17| ) one also obtains 



ln(I+ V-H" -1 ; 



(5.1) 



a D (H v ,H) = -K D ^ 2h (H) V 2 + 0(V 3 ) 



Now, taking the limit V 



K D - 2h (H) 



lim 

V->0 



0, one may show that 

res H- 2 , 



a D (H v ,H) 1 



4V 2 



Ah 



(5.2) 



(5.3) 



in agreement with Ref. [25|. 
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6 Some physical examples 



Here we are going to consider some basic examples in which the multiplicative anomaly may 
actually have physical consequences. Our aim has been, in fact, to look for the most simple 
situations (from the point of view of the number of space dimensions and of the order and 
nature of the differential operators involved) that might already exhibit a non-zero anomaly. It 
should not come as a surprise that, among them, the harmonic oscillator plays an important 
role. It yields non-trivial results easy to calculate analytically in any number of dimensions. 



6.1 Presence of the anomaly for Dirac-like operators in one space dimension 



Consider the square root of the harmonic oscillator obtained by Delbourgo in Ref. [26]. This 
example has potentially some interesting physical applications, for it is well known that a fermion 
in an external constant electromagnetic field has a similar spectrum (Landau spectrum) . Exactly 
in the same way as when going from the Klein-Gordon to the Dirac equation and at the same 
price of doubling the number of components (e.g., introducing spin), Delbourgo has constructed 
a model for which there exists a square root of its Hamiltonian, which is very close to the one 
for the harmonic oscillator. It is in fact different from the Dirac oscillator introduced by several 
other authors, corresponding to the minimal substitution p —> p — iar. The main difference lies 
in the introduction now of the parity operator, Q. Whereas creation and destruction operators 
for the harmonic oscillator, a = P ± iX, are non-hermitian, the combinations T>^ = P ± iQX 
are hermitian and 



(V 



±\2 



P 2 + X 2 tQ = 2tfosc T Q- 



(6.1) 



Notice that the parity term commutes with H sc- Doubling the components (cr, are the Pauli 
matrices) 



P 



-lO\ 



d_ 

dx ' 



X — ► a\x, Q — > cr 2 , 



the operators P ± are represented by 



V 



± 



• 9 + 

OX 



(6.2) 



(6.3) 



In the sequel, we will only consider the operator T> = T> + . It has for eigenfunctions and 
eigenvalues, respectively, 



-te 



-x 2 /2 



2 n+i ( n - 1)!0F 



H n .. 1 (x)±H n (x)/y/2ii 
H n ^(x)TH n (x)/V2^ 



X n = ±v2n, n > 1, 




V>oOr) = ^= , , A = 0, (6.4) 



where the H n (x) are Hermite polynomials. 

The two operators we shall consider for the calculation of the anomaly are T> and T)y = T>+V, 
V being a real, constant potential with \ V\ < y/2, that goes multiplied with the identity matrix in 
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the two (spinorial) dimensions (omitted here). Notice that T> and T>y are hermitian, commuting 
operators. The zeta function for the operator T> reads 



C(s\V) =Y,K S = Et 1 + (V^F = [1 + {-1)- s ]2- s, \r{s/2), (6.5) 



n=l 



Cr( s ) being the usual Riemann zeta function, which has a simple pole at s = 1. Furthermore, 



the manifold is not compact, thus, by direct use of Eq. ( 2.11 ) we readily get 



a(V,V v ) = ^--lnV. (6.6) 

The logarithmic term is due to the presence of a zero mode. 

How does this match with the Wodzicki formula? First of all, we point out that one has to 
deal with zero modes. Thus some care must be taken to properly treat them. Secondly, since 
we are working in a non-compact manifold, we cannot used the local formula to evaluated the 
non-commutative residue. A direct check, modulo the zero mode problem, that makes use of 
Eq. (|3.9|) , yields the multiplicative anomaly above. It thus seems that Wodzicki's expression 
requires only small modification in order to deal with spinorial operators as the ones we have 
here. 

We conclude by pointing out that we have here, before us, the first and most simple example 
of the presence of a non-trivial anomaly for operators of degree one in a space of dimension one 
(spinorial, however). 



6.2 Generalization of the Dirac-like operators to D dimensions 

Referring again to the work by Delbourgo in Ref. p6[ , the above operators T> + and V>~ admit 
possible extensions to D dimensions, which have for eigenvalues, respectively, 

A+ = ±y/2(2n r + l), X- = ±^2(2n r + l + D), n r , I = 0, 1, 2, . . . (6.7) 

Notice that T> + exhibits a zero mode (for n r = I = 0), what is not the case with D~ . Each of 
these two operators provides a different example for the calculation of the anomaly, and we will 
denote the corresponding partners by Vy (see the preceding subsection). 
The basic zeta functions are now 

oo ' 

C + (s\H) = [l + (-l)- s ] 2- s / 2 ^ {2n r + l)- s ' 2 

n r ,l=0 

= [i + (-!)-] 2-/ 2 C 2 (|,0| (2,1)), (6.8) 

oo 

C-(s\H) = [l + (-l)- s ] 2 - s / 2 ^ (2n r + l + Dy s/2 

n r ,l=0 

= [l + (-l)-] 2-'/%(i,D\ (2,1)). (6.9) 

As usually, the prime means in the first expression that the zero mode has to be excluded from 
the sum (i.e., the term with n r = I = 0). Here (,2{s, b\(a\, a 2 )) is the Barnes zeta function in 2 
dimensions. It has simple poles at the points s = 1, 2 (in general at s = N, N — 1, . . . for (jv) 
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and the residues are well known, as given by generalized Bernoulli polynomials at those points 
(see Refs. |27-29t|). Again using Eq. (|2.11|) , one gets 



a(V+,Vt 



a_ = a(T> , T>\ 



2_ — : Res C2 (j, d\ (2, 1)) 



i=i 



2i 



(6.10) 
(6.11) 



where, as before, the logarithmic term is due to the presence of the zero mode, and ip is the 
digamma function. For D = 2, the result is 



a+ 



a. 



y2 hiy2 
T \ 12 

4 1 12 



1 



1 



InV, 



Let us now check with the result obtained from Wodzicki's formula 

, 2" 



res 



(ln(^M 2 - ai ; 



2aia2(ai + 02) 



1 

- res 
4 



V 



(6.12) 



Here ai = 02 = 1. Looking for the non-commutative residue, making use again of Eq. (p.Sj), of 
course one obtains Eqs. Q6.10 ) and ( |6.11 ) (modulo the logarithmic term corresponding to the 
zero mode, which is immediate to supply). From the physical point of view, it is not completely 
clear if these operators make sense for any value of D. In fact, in Ref. [^] some doubts were 
arisen concerning their precise physical meaning in four and higher dimensions. 



6.3 Harmonic oscillators in D dimensions 

Let us recall the case of the harmonic oscillators in D dimensions, with angular frequencies 
(u>i, ... ,ojb)- The eigenvalues read 



Afj, = n ■ lo + b . 



n 



(ni,...,n D ), u = (ui, ...,u D ) , 



1 



D 



(6.13) 



fe=l 



and the related zeta function is the Barnes one Cd(s, b\ui), whose poles are to be found at the 
points s = k {k = D, D — 1, . . . , 1). Their corresponding residua can be expressed in terms of 



(6.14) 



generalized Bernoulli polynomials B D _ k (b\uj). They are defined by [p0| 



t D e -at 



njLiCi-e- 6 **) n i= i^ n=0 

The residua of the Barnes zeta function are then: 

\D+k 



Res Cd(^> b\w) 



(k-iy.(D-k)\UiLi^ 



k = D,D — 1,... 



(6.15) 
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Now, if V is a constant potential, from Eq. Q2.11| ) we easily obtain 
a(H,H v ) 



[ ^ h + 1>(D-2k)]B%\b\Q) y2k 



211?-!-. 



j=l"3 k=l 



{2k)\{D-2k)\ 



(6.16) 



Notice that in our case the generalized Bernoulli polynomials of odd order vanish: B[ D \b\ui) 
B ( 3 D \b\Lu) = ■■■ = 0, for any D. On the other hand, the remaining generalized Bernoulli 



polynomials are never zero, in fact 



B, 



(D), 



Bi D \b\u) 



D 



12 ^ 



i=l 



B[ D \b\Q) = 



1 

24 



D 



bI>?+E 



2 2 



1=1 



(6.17) 



B^\b\Q) 



5 

96 



31 x g 7 x 
70^^ + 7n ^ 



4 , , 2 



2 2 2 



i=l i<j<k 

As a consequence, the anomaly does not vanish in any case, not for D odd or D = 2, whatever 
the frequencies be. Moreover, only even powers of the potential V appear. Again, since the 
manifold is not compact, a direct use of Eq. (^) confirms the validity of the Wodzicki formula. 



7 Conclusion 

In this paper, we have considered the issue of the evaluation of functional determinants that 
typically are present in the one-loop approximation of the Euclidean partition function. We have 
stressed that the regularised determinant of a product of elliptic operators has a multiplicative 
anomaly factor. In general, such a factor is not vanishing and it should be taken into account in 
associated physical applications, since it can play an important role. The non-commutative case 
of the multiplicative anomaly formula has been discussed in some detail, making use of different 
techniques. We have started with an elementary application of standard perturbation theory of 
linear operators to a particular case involving a product of two of them. 

Then, more powerful techniques of ^DOs and the Wodizcki theory of the non-commutative 
residue have been employed in a useful re-derivation of the general multiplicative anomaly for- 
mula. 

One of the most interesting results obtained in the paper deals with the possibility to notably 
simplify the general formula for the non-commuting case in the lower dimensional situations: 
D = 2, D = 4, both of which turn out to be very relevant for physical applications. 

Several basic examples have been discussed with care, mainly in the non compact case. Here, 
another result of our paper is the formulation of a conjecture that extends the validity of the 
Wodzicki formula to the calculation of the multiplicative anomaly in the case of non-compact 
manifolds. 
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